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Abstract. We discuss the physics of the VtT- Limit, i.e. when the star is 
unbound as a result of both rotation and radiation pressure. We suggest that the 
f2r- Limit is what makes the Humphreys-Davidson Limit. Stellar filiations are 
discussed, with an emphasis on the final stages, in particular on the final masses, 
their angular momentum and the chemical yields with account of rotation. A 
possible relation between WO stars and GRB is emphasized. 



1. Introduction 

Massive stars are at the crossroads of major subjects in Astrophysics. They 
have a key influence on the spectral and chemical evolution of galaxies. Their 
contributions to the spectrum of starbursts is observed up to large redshifts. 
They are also the essential components of the metallicity Z = populations, 
the progenitors of WR stars, supernovae and gamma ray bursts (GRB). 



2. The or Limit: physics of the interaction of mass loss and rotation 

The physics and evolution of massive stars is dominated by mass loss and by 
rotational mixing. At the origin of these two effects, we find the large ratio 
T/p of temperature to density in massive stars. This enhances the ratio of 
radiation to gas pressure, which goes like T^/p and favours stellar winds. A 
high T/p enhances rotational mixing, since mixing by shear turbulence scales 
as the thermal diffusivity K = AacT^ /3cp k/?^. Also, the velocity of meridional 
circulation scales as the ratio L/M of the luminosity to mass. 

Let us recall that there are 3 main effects of rotation. 1) Structural ef- 
fects due to the centrifugal force: these are small in the interior, but may 
produce a large distorsion at the surface. 2) Rotational mixing: this pro- 
duces internal transports of chemical elements and angular moment um. Mixing 
by shears is most efficie nt for the transport of chemical elements ( Zahn 19921 : 



iMevnet Sz Maeder 200lh . while meridional circulation dominates the transport 
of angular momentum. 3) Enhancement of mass loss by rotation: an impor- 
tant case is the fiP- Limit, where the star reaches break-up as a result of both 
high radiation pressure and rotation. The von Zeipel theorem states that the 
local radiative flux F of a rotating star is proportional to the effective gravity 
Qes- Th us, there is a much larger flux and a higher Tefj at the pole than at the 
equator ((Mevnet Sz Maeder 19971 ). This latitudinal dependence of Teff leads to 
asymmetric mass loss and also to enhanced average mass loss rates. 
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Figure 1. Left: Iso-mass loss distribution for a 120 Mq star with 
log = 6.0 and T^g = 30000 K rotating at a fraction 0.80 of break-up 
velocity. Right: The same with Toff = 25000 K. 



On a rotating star, one must consider the flux F{'d) at a given colatitude 
'& as given by von Zeipel's theorem, F(i!}) = ~ 4^ga/ geff[l + C(^)]- The term 
C{'&) is in general negligible. The Eddington factor is a local quantity rQ(i?) 
depending on and rotation. We define it as the ratio of the local flux F{t!)) 
given by the von Zeipel theorem to the maximum possible local flux, which is 
-^iim('!9) = --^^Oesi^)- Thus, One has 

FuU^) A^cGM (l - ^) 

where the opacity k{'&) depends on the colatitude ■(?, since T^g itself depends on 
i9. For electron scattering, k is constant and rf2(i?) is the same at all latitudes, 

i.e. = r/ ^1 — 27rG'pm ) ' where F is the usual expression. Eq. (1) shows that 
the maximum luminosity of a rotating star is reduced by rotation. We see that 
the dependences of F{'&) and of Fii^{i!}) with respect to g^g have cancelled each 
other in the expression of Fq(??). Thus, if the limit Ff7('!9) = 1 happens to be met 
at the equator, it is not because ^efF is the lowest there, but because the opacity 
is the highest! The theory of radiative winds applied to a rotating star leads to 
an expression of the mass flux as a function of colatitude ( Maeder 1999t ). Figs. 1 



left and right (Maeder & Desiacqucs 2000) illustrate the distribution of the mass 
loss rates around a massive star of 120 Mq for two different Tefj. For a star hot 
enough to have electron scattering opacity from pole to equator, the iso-mass 
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Figure 2. Left: simulation of a short shell ejection by a massive star with 
anisotropic mass loss as in Fig. 1 left. Right: the nebulae around AG Carinae 
l|Nota fc Clampin 1997j) . 



loss curve has a peanut-like shape (Fig. 1 left). This results from the fact that 
the pole is hotter ( "g'eff-effect" ) . For a rotating star with a lo wer Tpff (Fig. 1 
right), a bistability limit i.e. a steep increase of the opacity LLamers et al. 1995^ 
may occur somewhere between the pole and the equator. This "opacity-effect" 
produces an equatorial enhancement of the mass loss (Fig. 1 right). In Fig. 2 left, 
we show the model of a short shell ejection with mass l oss corresponding to the 
peanut-shape. The corresponding image of AG Carinae (|Nota fc Clampin 1997^ 
is shown in Fig. 2 right. 

The anisotropics of mass loss influence the loss of angular momentum. Polar 
mass loss removes mass but relatively little angular momentum. This has a 
great incidenc e on the evolution of the most massive stars with high rotation 
( Maeder 200^ . The mass loss rate M(Q,) of a rotating star compared to that 
M(0) of a non-rotating s tar at the same location in the HR diagram is given by 
( Maeder fc Mevnet 200(1 . 
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where a is a force multiplier ( Lamers et al. 19951 ). For a 10 Mq star on the MS, 

may reach 1.5. For the most luminous stars which have a value F close 

to 1.0, this may be orders of magnitude, when the star reaches break-up at the 
OF-Limit. 

Often, the critical velocity in a rotating star is written as = ^^(1 — F). 
This expression is not correct, as it applies only to uniformly bright stars. In- 
deed, the critical velocity of a rotating star is given by the zero of the equation 
expressing the total gra yitv qt7,t, = q^av + qrnt. + ftad = g7s [1 - (■!?)]• This 
equation has two roots ( Maeder Mevnet 20001 ). The first that is met deter- 
mines the critical velocity. The first root is as usual fcrit,i = (f^^)^' where 
Rph is the polar radius at break-up. The second root ?;crit,2 applies to Eddington 
factors bigger than 0.639. It is equal to 0.85, 0.69, 0.48, 0.35, 0.22, times fcrit,i 
for F = 0.70, 0.80, 0.90, 0.95, 0.98 and 1.00 respectively. 



3. Close to the Eddington Limit, evolution unavoidably reaches the 
rJF Limit 

On the whole, we may distinguish 3 critical cases, where the outer layers escape: 
1.- The usual F-Limit, when radiation effects largely dominate over rotation; 
this is the classical case. 2.- The il-Limit, when rotation effects, rather than 
radiative effects, are determining break-up. This is relevant for Be-stars, which 
are far enough from F equal to unity. 3.- The fiF-Limit, when both rotation 
and radiation are important. This is the case which applies to the most massive 
stars. Even for a rather small initial rotation velocity of, say, > 50 km/s, a star 
with a high F will reach critical velocity during its MS evolution. This occurs 
because the rotation velocity increases during MS evolution (see Fig. 3), while 
the critical velocity, which is given by Ucrit,2) decreases. The mass loss rates 
increase strongly (cf. Eq. 2), when the cr itical OF-Lim it is approached. Such 
a situation has already been considered ( Langer 1997t ). however for the case 



of solid body rotation, while here the evolution of the angular momentum is 
followed with appropriate critical velocity and mass loss expressions. 

The evolution of star reaching the J7F-Limit is illustrated in Fig. 3. On 
the left, it shows different cases of tracks in the HR diagram. The track with 
■Willi = 500 km/s and an average mass loss rate M of 10~^-^ Mq -yr'^ is calculated 
according to r ecent M-rates (|VinketaOoo3) . With these rates, even the zero 



rotation track would turn bluewards! Thus, at I2OM0, these rates are likely 
too high (at least for producing a star like rj Carinae). With rates a factor of 
2-3 lower, a bluewards evolution is obtained only for a very fast rotation, such 
as 650 km/s. For initial rotation Vini = 600 km/s or lower, redward evolution 
is obtained, with tracks having a different MS extension (for low rotation, the 
tracks are higher and extend further in the HR diagram than the track for zero 
rotation; for fast rotation, the extension is shorter as the result of the lower 
opacity due to higher helium enrichment). All tracks with Vini > 50 km/s reach 
the QF-Limit with break-up velocity, thus experiencing phases of extreme mass 
loss. 
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Figure 3. Left: Evolution of a 120 Mq star with different mass loss rates and 
initial rotational velocities Uini, see text. Right: Evolution of the rotational 
velocities of models with different Vini- The envelope of the various curves is 
the location where the stars reach their critical velocities. The models with a 
lower Vini reach the critical velocity later. A model with solid body rotation 
with Willi = 200 km/s is also shown. It reaches the critical velocity faster. 

Fig. 3 right illustrates the corresponding behaviour of the rotational ve- 
locities. We see that the velocities increase during evolution. The faster 'I'jni? 
the earher the Qr-Limit is reached (envelope of the curves). We see that solid 
body rotation, which is the extreme case of coupling of the angular momentum, 
leads the star very promptly to the OF-Limit. The milder coupling realized by 
meridional circulation leads the star to the JIF-Limit more slowly. We conclude 
that the most massive stars reach an r^F-Limit, which is different for each star 
according to its rotational velocity. The fiF-Limit is the physical reason for the 
observed Humphreys-Davidson Limit (,HumDhrevs &: Davidson 1979. ). 

4. What occurs at the fiF Limit ? 

For now, hydrodynamical models of LBV outbursts do not exist, nevertheless 
we can make a few remarks on their properties: 

-1. Likely, the most massive stars reach the r2F-Limit on their redwards track 
as illustrated in Fig. 3 left. The slightly less massive stars, firstly evolving to the 
red, may likely reach the Limit on their bluewards tracks, because as a result of 
their mass loss, their L/M is larger when they go back to the blue. Also, due to 
the coupling of angular momentum by convection, they may rotate faster when 
evolving bluewards (Langcr 1998). 

-2. As stated above, the M-rates increase as given by Eq. (2), when the QT- 
Limit is approached. However, at the Limit, the M-rates are no longer given by 
the stellar wind theory. We think that the M-rates are such that the star just 
remains at the critical velocity, (or in practice, at a fraction e.g. 0.98 or 0.99 of 
it). If the M-rates would be larger the star would soon become subcritical, if 
they would be smaller, the star would rotate much above critical. 
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Figure 4. Left: The broken line shows the V-Ught curve for a star with 
logTeff = 3.8 experiencing a mass ejection of I.OMq, the continuous hne 
shows the same for an ejection of 3 Mq , the last part of this curve would be 
similar to the broken line. Right: Evolution of the massive stars in the plot 
of the surface H-content Xs vs. luminosity. 



-3. Some layers become supra-Eddington, which produces a density inversion 
and convective instability ap pears. Thus, even a B-type supergiant can have 
convection in its outer layers ( Maeder 198Ci ). 

-4. In the very massive stars, the thermal timescale may become smaller than the 
dynamical timescale in the outer layers. This means that the thermal structure 
of these layers can readjust during a dynamical event. Thus, the ionization 
front may move inwards and allow a substantial mass fraction to participate 
to a dynamical event. This i s the idea of the "geyser model" (fMaeder 19921 : 
iHumphrevs fc Davidson 19941 : iMaeder 1997l 'l. which leads to an estimate of the 
amount of mass ejected as a function of the lum inosity. This has be en confirmed 
by IRAS observations of dust in LBV nebulae ( Hutsemekers 1997h . 
-5. As a result of the mass ejected in the outburst, the star is shifted to the blue 
in the HR diagram, then it slowly recovers its previous Tefj . The luminosity keeps 
nearly constant, apart from the energy lost in the mechanical event. Simulations 
of sudden large mass ejection show several properties: - a) The extension of blue 
shift in the HR diagram depends on the ejected mass AM; - b) The timescale 
of the blueshift is determined by the mass loss rate in the outburst; - c) The 
light curve in V magnitude (cf. Fig. 4 left) is mainly determined by the change 
of the bolometric correction (see lightcurve of rj Car); - d) The recovery time, 
i.e. the time for the star to reach its former equilibrium depends on the ejected 
mass AM. It is about 125, 350 and 750 yr for AM = 0.3, 1.0 and 3.0 Mq. 



5. Filiations and further evolution of the massive stars 

There are various evolutionary sequences: for the highest masses with mass loss, 
the stars may keep to the blue. Then in an interval of lower masses, the tracks 
have a certain extension to the red and even in a lower mass interval the full 
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Figure 5. HR diagram of massive stars with metallicity Z=0.02 



extension to the red supergiants is present (Fig. 5 ). We may distinguish the 
following tentative filiations at solar Z: 

• Of - WNL - (WNE) - WCL - WCE - SN (Hypernova low Z ?) 

- Of/WNL< - >LBV - WNL(H poor)- WCL-E - SN(SNIIn?) 

- BSC - LBV < - > WNL -(WNE) - WCL-E - SN(SNIb) 
- WCL-E - WO SN (SNIc) 

- BSC - RSG - WNE - WCE - SN(SNIb) 
OH/IR < - > LBV ? 

-(BSC)- RSG - BSC (blue loop) - RSG - SN(SNIIb, SNIIL) 

• RSG - (Cepheid loop, M < 15 Mq) RSG - SN (SNIIL, SNIIp) 

> means back and forth motions between the two stages. The 
limits between the various scenarios depend on metallicity Z and rotation. The 
various types of supenovae are tentatively indicated. 

The variation of the H-surface content as a function of luminosity is a 
constraining test (Fig. 4). The most massive stars go almost vertically down . 
The LBV stars, according to a recent collection of data ( Stothers &: Chin 20001 ) . 
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are located between log L/L© = 5.4 and 6.2 with Xg between 0.30 and 0.40. The 
WN stars are in almost the same range of L while their Xg are between 0.45 and 
0. This overlap supports the back and forth transition between the WNL stars 
and the LBV. The further evolution from LBV leads them to WN stars, firstly 
WNL and then WNE. The W R lifetimes as a function of initial mass, rotation 
and Z have been calculated ( Mevnet Sz Maeder 20041 ). The comparison with 



observations is considerably improved when rotation is accounted for. 

Fig. 6 left shows the relative enrichments in N/C for stars of 9 and 40 Mq 
at different Z. In the lower mass domain, the enrichments are higher for stars 
with lower Z since mixing is in general stronger. This results from the steeper 
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Figure 6. Left: Relative enrichments in N/C for rotating stars of 9 and 40 
M© at different Z but with the same Wini ~ 300 km/s. Right:The final masses 
at the time of SN explosion as a function of mass and Z . 



internal fi-gradients in lower Z stars ( Mevnet &: Maeder 200lh . For masses 
above 30 Mq, mass loss dominates the evolution and as the M-rates are larger 
at higher Z, the enrichments are also larger at higher Z . 



6. Final stages in relation with 7-ray bursts (GRBs) 

As a result of mass loss, the stellar masses decrease very much. Fig. 6 right 
shows the final masses at the time of supernova explosion as a function of the 
initial masses and Z. For solar Z or higher, all stars with initial masses above 
20 Mq finish their life with masses of about 10 Mq. At lower Z, the final masses 
are higher. It should not be concluded that for Z = the final masses are equal 
to the initial ones, because such stars are likely to spend a fraction of their MS 
lifetime at break-up and loose mass (see Meynet, this conference). The different 
final masses lead to different types of supernovae. 

Fig. 7 left shows the distribution of the specific angular momentum j at 
various stages of an initial 25 Mq star. The final distribution of j is essen- 
tially shaped during MS evolution. Thus, the treatment of the transport of 
angular momentum in the MS phase i s an essential aspect . The various models 
(|Heger et al. 20001 : iHeger et al. 2003 I : iHirschi et al. 20041 ) lead to rather simi- 
lar distributions, despite their differences in input physics. A specific angular 
momentum above 10^^ cm^ s~^ is ne cessary for the collapsar model to work 
( Heger et al. 2000l : lHeger et al. 200'3Tl and account for the (GRB). These mod- 
els fulfill this condition, however a lot of angular momentum has to be lost in 
the explosions in order to account for the rotation rate of pulsars. 

The association of GRB with hypernovae of the class of SNIc is sup ported 
by several observations f Mazzali et al. 200^ : iPodsiadlowski et al. 200 J ). SNIc 
result from the explosion of a star without H and with little or no He. This 
corresponds to a rare category of WR stars: the so-called WO stars. WO 
stars ares subsequent to WC3 objects, they show the products of He-burning, 
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Figure 7. Left: Distribution of the specific angular momentum in a model 
with an initial mas s of 25 Mm. The distribution is shown at various evolution- 
ary stages ( Hirschi ~t al. 200^ . The distribution in the last model coincid e 
well with results from Heger et al. l|Heger et al. 200(7: ' Heger et a l. 200.3^. 
Right: Evolution of the ratios (C-|-0)/He as a function of the luminosity 
at the surface of 60 and 120 Mq rotating models for various initial metal- 
licities (see text). Long-short dashed curves show the evolution of Z = 
0.004 models, dashed-dotted curves, short-dashed curves and continuous lines 
show the evolutions for Z = 0.008, 0.020 and 0.040 respectively. The cor- 
respondence between th e (C-|-0)/He ratios and the different WC subtypes 
ijSmith fc Maeder 1991(1 is indicated on the right. 




Figure 8. Left: the yields x IMF for models without rotation. Right: the 
yields x IMF for models with Vmi = 300 km/s l|Hirschi et al. 2004) . 



with an excess of C+O with respect to He and O > C. WO stars are rare 
and result from t he evo lution of stars with M > 60 Mq at low metallicity only 
ISmith & Macdc r 199ll ). These authors show that early WC types and WO 
stars are found in lowe r Z regions. This may seem surprising. The reason is 
(|Smith fc Maeder 199ll 'l: at high Z, mass loss is high, thus when the products 
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of the 3a reaction appear at the surface, they are in an early stage of nuclear 
processing, i.e. with a low C + O /He ratio (see Fig. 7 right). This corresponds 
to late WC stars. At low Z, the products of the 3a reaction rarely appear at 
the stellar surface and if they do it (e.g. in case of high rotation) this occurs 
very late in the evoluti on, i.e. when C + O/He is high and this is a WO star 
( Smith &: Maeder IQQlI 'l. Thus, the rare WO stars may be the progenitors of 
SNIc and GRB. 

Nucleosynthesis is also influenced by rotation ( Hirschi et al. 20041 ). Below 
30 Mq, the a-elements are enhanced due to the larger cores. Above 30 M©, 
mass loss is the dominant effect and more He is ejected before being further 
processed. Fig. 8 shows the yields multiplied by the initial mass function. Due 
to the weighting by the IMF, the production of oxygen and of a-elements is 
globally enhanced, while the effect on the He-production in massive stars is 
rather small. 
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